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Based on recent magnetic-quantum-oscillation, ARPES, neutron-scattering and resistivity data, 
we propose that the extraordinarily high superconducting transition temperatures in the cuprates 
are driven by an exact mapping of the incommensurate spin fluctuations onto the d^i_y-2 Cooper- 
pair wavefunction. Strong evidence for this comes from the direct correspondence between the 
incommensurability factor 5 seen in inelastic neutron scattering and the inverse superconducting 
coherence length (see Fig. 2). Based on these flndings, one can specify the optimal characteristics 
of a solid that will exhibit "high Tc" superconductivity. 



In spite of strong evidence for an order parameter ex- 
hibiting d^2_y2 symmetry [l| in the "High Tc" cuprate 
superconductors, experiments establishing a direct link 
with the pairing mechanism have not been forthcom- 
ing. However, with the advent of magnetic-quantum- 
oscillation measurements on underdoped cuprates 
0, 01 and consequent proposals for the Fermi-surface 
topology that are consistent with experiments ranging 
from neutron scattering to ARPES [1, [1, 0], there is 
the opportunity to resolve this issue. We find that 
the dx2_y2 Cooper-pair wavefunction is exactly tuned 
to the real-space form of the fluctuations responsible 
for the incommensurate antiferromagnetic scattering ob- 
served in inelastic neutron experiments on all doped 
cuprates 0, i, B K Ell, El (see e.g. Fig. [T]i). The 
dispersion of these fluctuations, which are likely to be 
driven by the varying topology of the Fermi surface as 
the hole density p increases [J] , provides a very plausible 
energy scale for Tc in the cuprates [1]. 

Most treatments of superconductivity consider the 
state as a condensate 13, 13, El, EE E3|; such a pic- 
ture, it is natural to treat the Cooper-pair wavefunction 
in fc-space. However, we are interested in the similar- 
ity of the spatial distribution of spin-fluctuations and 
the Cooper pair, and so we consider the form of the 
two-dimensional d^^-y^ Cooper-pair wavefunction in real 
space }ll.ll7|: 



V;(r) cx cos(rfcF)(a;^ - J/^)e 



(1) 



Here r — ^JsP+y^ is the cylindrical polar radius, x 
and y are corresponding Cartesian coordinates, ftp is 
the Fermi wavevector 0, Q and is the supercon- 
ducting coherence length. Figure mb) shows a contour 
plot of ■(/'(r); the known diagonal nodal regions P, EB| 
are clearly visible. In such a plot, the lengthscale over 
which V'(r) is non- negligible is defined by (Eq- E])- 
We now consider similarities between this wavefunction 
and the real-space topology of the incommensurate an- 
tiferromagnetic fluctuations seen in inelastic neutron- 
scattering 0, i, i, E3, El El • 

Low-energy inelastic neutron-scattering data for both 
La2_^Sr^Cu04 SEl and YBaaCugOr-^ [3 EH follow 
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FIG. 1: (a) Ener gy ( E) versus for several YBa2Cu307_a; 
compositions p^ . |TJ , with fits to Eq. [5] to obtain extrapo- 
lated values oi 5 &t E = 2 meV (dotted line), (b) Contour 
plot of d-wave Cooper-pair wavefunction in two-dimensions 
(Eq. [TJ; we assume a Fermi wavevector fcp = 7r/\/2a where 
a is the lattice parameter, while p = 0.1 for which we select 
^0 = 27 A from Table|l] (c) Simulation of the incommensurate 
neutron scattering peaks obtained by Fourier transforming 
Eq. m 5 and 5 correspond to p = 0.1 (Table ID . (d) Real- 
space spin-fluctuation map corresponding to the set of four 
incommensurate peaks in (b), also for p = 0.1 (see Table|T]|. 
We consider the simple case where sq is the same for all Q. 



an approximately inverted parabolic form 



a^(q - Qo 



i?(q) = £^0 1 - 



(2) 



where a is the in-plane latticeparameter and Eq is a 
doping-dependent energy scale 0, H, 0, E3, El ; some ex- 
amples are shown in Fig. [ija). Here, 5 depends on the 
composition and doping of the cuprate involved H, E3] ; 
some values are given in TableU ksE ^0, the brightest 
scattering intensity typically appears at a cluster of four 
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incommensurate peaks at Q = (±tt / a, zt(l±2d)TT / a) and 
(±(1 ± 25)7r/a,±7r/a) 0, i, B 0, [III, 0] ; a simulation 
is shown in the positive , positive qy quadrant of two- 
dimensional fc-space in Fig. [Ijc). The form of the scat- 
tering peaks in Fig. [T]^c) is already very suggestive of the 
dr^2_y2 wavefunction in Fig.^b); to make a quantitative 
comparison, however, we require a real-space representa- 
tion of the corresponding spin fluctuations. 
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FIG. 2: a The approximately linear correspondence between 
d = a/25 and for La2-a;Sri:Cu04 (o) and YBa2Cu307-a; 
(□). The line represents the ratio £,o/d — 3/2. b Plot showing 
S and £,Q^ in La2-a^Sri,Cu04 exhibiting a maximum versus p 
in a similar fashion to Tc (scaled for comparison), though 
falls more rapidly with p beyond optimal doping. 

A simple model of spin fluctuations that can produce 
the observed incommensurate scattering peaks is a sinu- 
soidal variation of the staggered moment modulated by 
an exponential damping factor; the latter term represents 
the fact that the antiferromagnetic fluctuations possess a 
finite correlation length ^ 7, 18]. The spatially- varying 
moment is thus 



sir,t) 



Q 



SQexp(-- +icjt)cos(Q.(r-ro)), (3) 



where lo is the angular frequency of the fluctuations, 
ro — (±rf, 0) or (0, id) with d — a/26, and the sum 
in Q runs over the values given above. The simulation 
shown in Fig. [TJc) is a Fourier transform of Eq. [3] that 
reproduces the gene ral form of the experimental neutron 
data BBBEinilil very well. Note that the choice of 
ro is quite critical; any other value would give significant 
intensity at q = Qo, at variance with the experimental 
spectra. This is an important point, to which we will 
return. Finally, Fig. [IJd) shows a time-averaged contour 
plot of Eq. [31 The similarity between the spin fluctua- 
tion distribution and the dx-2_y2 wavefunction (Fig.[ljb)) 
is most marked: both show similar angular and radial 
distributions. Less obvious, but equally germane, is the 
phase of each function plotted. The choice of ro neces- 
sary to reproduce the neutron data means that there is a 
TT difference of phase between adjacent lobes of the spin- 
fluctuation distribution, just as there is a tt difference in 
phase between adjacent lobes of the dr^-i_y-2 Cooper-pair 
wavefunction [ll. Il6l|. 

As mentioned above, the lengthscale of the Cooper- 
pair wavefunction is ^o! the corresponding lengthscale 




FIG 



3: Comparisons of the spin density amplitude and 
y2 wavefunction at selected values of p. The first col- 



umn shows plots of the incommensurate diffraction peaks ob- 
tained on re- Fourier transforming s(r) calculated using Eq. [3] 
and the published values of 5 and ^ listed in Table IH The sec- 
ond column shows the calculated |s(r)| while the third column 
shows Pc(r) = |*I'(r)p calculated according to the values of 
listed in Table |T] Rows 1, 2 and 3 consecutively correspond 
to p « 0.06, 0.1 and 0.152. 



for the spin fluctuation distribution (Eq. [3|) is d = 
a/25 0, S S [IS IHI- To further establish the causal 
relationship between the spin fluctuations and the super- 
conducting pairing mechanism, we now compare indepen- 
dent estimates of these lengthscales for several cuprates. 
The values of d are obtained from inelastic neutron- 
scattering experiments. In the La2_2;Caa;Cu04 cuprates, 
the dispersion relationships (Fig. [Ijb)) have been mea- 
sured down to low energies, providing accurate values 
of (5 0]. However, there is a loss of intensity at lower 
energy transfers in YBa2 Cu.'^ 07-^i possibly because of 
their greater homogeneity [lO, [lH , necessitating a down- 
ward extrapolation of the dispersion curves to obtain 5 
in the limit i? — > 0. The coherence lengths are in- 
dependently derived from magnetoresistance and other 

3, 



data [191 . [20| ; occasionally some interpolation is required 
to obtain values for the same hole doping and compo- 
sition as the neutron-scattering data. Fig. shows the 
corresponding values of i^o plotted against d. Table [J lists 
the values (plotted against doping for La2-a;Sr2;Cu04 in 
Fig. [lb). 

Remarkably, Fig. [2^ shows that the i^o versus d values 
for all compositions and dopings listed in Table [Tj lie on 
a single straight line £^o/d — 3/2, irrespective of whether 
La2_3;Sr3;Cu04 or YBa2Cu307_a; is considered, suggest- 
ing a common origin for the spin fluctuations and the 
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superconductivity in all of the cuprates. 

Having established the experimental relationship be- 
tween ^0 and d, we can now plot the evolution with 
increasing hole density p of the spatial distribution of 
the spin fluctuations (Eq. [3]) alongside the correspond- 
ing Cooper-pair wavefunction (Eq. [T]) using experimental 
values of both parameters; Fig. [3] shows the result. At 
small values of p, the incommensurate neutron-scattering 
peaks occupy a small area of fc-space SESIlli; the cor- 
responding real-space spin fluctuation distribution occu- 
pies a large does the Cooper-pair wavefunction. 
As p increases, the incommensurate peaks spread out in 
K-space @, Efl, El; hence, the spin-fluctuation spatial 
distribution is compressed, as is the Cooper-pair wave- 
function. 



TABLE I: Parameters for various cuprate superconductors, 
including hole doping p. For La2-2,Sr2,Cu04, S values mea- 
sured at neutron transfer energies of _B ~ 2-3 meV are taken 
from Refs. [3, [13]. For YBa2Cu307-a;, S values at the equiv- 
alent energy are taken from plots such as Fig. [2^ using 
data from Refs. [lol . [Til ]. Correlation lengths C are taken 
from Ref. M] while BCS coherence lengths Co are taken from 
Refs. (19I, [20(1 ■ occasionally requiring interpolation between ad- 
jacent compositions for precise doping matches to S. 
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Before turning to the origin of the incommensurate 
scattering peaks in the neturon data [IiHjB&Ej 11.1^. 
it is worth summarizing two factors that lead us to pro- 
pose a causal relationship (or one might say spatial reso- 
nance) between superconductivity and antiferromagnetic 
fluctuations in the cuprates. 

First: on setting Tq — (id, 0) or (0, ±d) where d = a/25, 
the TT difference in phase between the adjacent 'spin clus- 
ters' in Fig. [TJi is aligned with the it difference in phase 
between adjacent lobes of the dx2-y2 Cooper-pair wave- 
function. The spin fluctuations are therefore "in phase" 
with the Cooper-pair wavefunction. 
Second: independent of our choice of wavefunction, there 



exists a very direct correspondence between which de- 
termines the overall size of the Cooper pair, and 6, which 
quantifies the physical separation d = a/26 between ad- 
jacent clusters of spins in real space (see Fig. [2]). On 
combining both of these observations, the £,o/d — 3/2 
slope in Fig. [2^ can be seen to originate from a direct 
spatial overlap between the vector locations (±d, 0) and 
(0, ±d) at which the spin fluctuation amplitude is largest 
and those (±|^o, 0) and (0, ±|Co) at which the probabil- 
ity amplitude of the dx2_y2 wavefunction given by Eqn. 
([T]) is largest. While there are certainly systematic and 
random errors in the methods used to extract ^0 val- 
ues [1^ [i^l, they are consistent in La2-a;Sra;Cu04 with 
estimates obtained from strong-coupling variants of the 
BCS formulae Co = ^wf/ttDo and 2no/fcBTc « 4 - 5, 
where Dq is the T = order parameter 

Bin, 

usmg 

typical Fermi velocities vp « 6-8 x 10 ms""'^ determined 
from quantum oscillation experiments on YBa2Cu408 
and YBa2Cu306.5 0, This provides a natural ex- 
planation for the observed linear dependence of S on 
Tc A slightly stronger couphng^^ //cbTc 6 is 



required for the YBa2Cu307 

Incommensurate diffraction peaks in metallic systems 
are almost always associated with Fermi-surface nest- 
ing, in which the periodicity of spin or change modu- 
lation is determined by the topology of the Fermi sur- 
face [2l|. Whilst the tight-binding calculations of the 
cuprate Fermi surface usually give a single, large hole 
pocket d, 0], the recent experimental observations of 
(multiple) small pockets in the underdoped cuprates 0, 
[1] supports the suggestion that some form of nesting 
occurs 0, S 01- Added weight is given by the similar- 
ity of the incommensurate mode dispersion (Fig. [Hb)) 
to those observed in itinerant antiferromagnets such as 
Cr and V2 



.O3 



23 



24j (both above and below the 
Neel temperature). In such a picture, the cuprate Fermi 
surface plays an increasingly important role as the hole 
doping p is increased. Thus, just as in V2-a;03 
large-moment antiferromagnetic insulator behavior in the 
cuprates eventually gives way to small moment incom- 
mensurate itinerant antiferromagnetic behavior as the 
system becomes more metallic [a, Consistent with 
itinerant magnetism, the orbitally-averaged Fermi veloc- 
ity tip = y/2ehF/m* « 8 x 10 ms~^, where F is the 
Shubnikov-de Haas oscillation frequency, of the pockets 
in YBa2Cu306.5 3 comparable to the mode velocity 
vo = 2aE,/'Kh5 « 14 X lO^* ms'^ (at E = Q) that one 
obtains fitting Eqn. [5] to the i?-versus-q data points ob- 
tained from inelastic neutron scattering experiments on 

(Fig. mh)) 



11, m 



As 



the same sample composition 
p increases toward optimum doping, Eq in Fig. [TJa) also 
increases 0, H , providing a suitable characteristic en- 
ergy scale 10s meV, i.e. ~ 100s of Kelvin) for Tc 3]. 

In conventional BCS superconductors, the quasipar- 
ticle interactions that result in pairing are via charge 
coupling to acoustic phonon modes [ij, [l5| . The in- 
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commensurate spin fluctuations have a dispersion rela- 
tionship (Fig. [TJa)) analogous to the acoustic phonons 
{i.e. approximately linear as E —^ 0, saturating at a 
maximum energy that defines the energy scale of Tf.). 
The remaining ingredient in the problem is therefore a 
coupling mechanism between the spin fluctuations and 
the charge inherent in the Cooper pair. This is found in 
the large onsite correlation ener gy U , which inhibits dou- 
ble occupancy of spins or holes [26| . Consequently, local 
variations in the density of holes Aph and spin-density 
amplitude As are expected to be subject to the relation 



Aph oc — A|s| 



(4) 



Pairing mechanisms involving this behavior have been 
considered both in the weak (small Hubbard U) 
and strong (large Hubbard U) [28| coupling limits, al- 
though typically in conjunction with long-range antifer- 
romagnetic order. However, there is no reason to sus- 
pect that such mechanisms will not apply in regimes 
where the antiferromagnetism is strongly fluctuating [t^ . 
Therefore, because of the reciprocity relationship (Eq.[4]), 
the slowly varying modulation of the spin fluctuation in- 
tensity should be accompanied by a concomitant charge 
modulation /Oh(r) oc — |s(r)|. The very similar form 
of Ph(r) oc — |s(r)| in the second column of Fig. [3] to 
Pc (X |^'(r)p in the third column of Fig. [3] therefore pro- 
vides direct evidence for a 'spatial charge commensura- 
bility' between the Cooper pair wavefunction and incom- 
mensurate spin fluctuations. 

Qualitatively, the Cooper pairs in this work are spa- 
tially compact near optimum doping (Fig. [3]), resembling 
the stron g co upling spin-bipolarons of Mott and Alexan- 
drov 2^ 3]- Away from optimum doping, they become 
spatially extended, like a weak-coupling d^2_y2 variant 
of the 'spin bags' proposed by Schrieffer et al. [23|. The 
difference between these pictures is that the fluctuations 
mediating the superconductivity in the present proposal 
are characteristic of a system on the brink of long-range 
order, rather than an established antiferromagnet. 

In summary, based on measurements including Fermi- 
surface studies neutron scattering data (tI. Isl. [glliol. 



191 120] , we propose that [19] 



11| and resistivity experiments 
the unusually high superconducting transitions in the 
cuprates are driven by an exact mapping of the incom- 
mensurate spin fluctuations onto the d^2_y2 Cooper-pair 
wavefunction. The spin fluctuations are driven by the 
Fermi-surface topology, which is prone to nesting [3, @] ; 
they couple to the itinerant holes via the strong on-site 
Coulomb correlation energy U, which prohibits double 
occupancy of spins or holes [26| . The maximum energy 
of the fluctuations (~ 100s of Kelvin 0, H, El) gives 
an appropriate energy scale for the superconducting T^. 
Based on these findings, one can specify the features nec- 
essary for a solid to exhibit "high Tc" superconductivity; 
(i) the material should be quasi-two-dimensional, with 
the conducting layers exhibiting four-fold symmetry (to 



ensure that fluctuations are optimally configured to a d- 
wave order parameter); (ii) the material should have a 
Fermi-surface topology susceptible to nesting, so as to 
produce antiferromagnetic fluctuations with a large de- 
gree of incommensurability S] (iii) however, the electron- 
phonon coupling should be moderate, to prevent forma- 
tion of stripe or charge-density- wave- like phenomena that 
compete with supercoductivity [2^; in other oxides, the 
larger electron-phonon coupling dominates, preventing 
superconductivity [30| . 
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